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1. Introduction 

A classic puzzle asks for the resistance between vertices at the ends of a long diagonal when the edges of 
a cube are replaced by 1-ohm resistors. The solution relies on the observation that for each of the endpoints, 
the three adjacent vertices are at the same potential, by the symmetry of the cube under a 120° rotation 
about the long diagonal. The network is thus equivalent to one in which three resistors in parallel are in series 
with six resistors in parallel and with three resistors in parallel, for a total resistance of 1/3 + 1/6 + 1/3 = 5/6 
ohms. (This problem seems first to have appeared in 1914 in a book by Brooks and Poyser [B4].) A natural 
question is: what happens when the 3-dimensional cube is replaced by an n-dimensional hypercube? 

Reasoning as before, we observe that all the vertices at a given distance from one of the endpoints of the 
long diagonal are again at the same potential, so the network is equivalent to a series connection of parallel 
connections of resistors. Since there are n + 1 distances (0,1,..., n) from one endpoint, there are n parallel 
connections. There are (^) vertices at distance k from the endpoint, and for < fc < n — 1, each of these 
vertices is connected by n — fc resistors to vertices at distance fc + 1 . Thus the total resistance is 

E 



0<fc< 



n 

Q<k<n 



where we have used the identity {^)in ^ k) = n("j^^), which is easily seen using the expressions for binomial 
coefficients in terms of factorials. (This n-dimensional version of the problem was posed in 1976 by MuUin 
and Zave as Problem E 2620 in the American Mathematical Monthly [M], with a solution by Jagers [J], and 
again in 1979 by Singmaster as Problem 79-16 in SIAM Review [SI], with a solution by Rennie [R].) 
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The numbers being summed in (1.1) are elements of the "harmonic triangle", considered by Leibniz as 
a companion to the "arithmetic triangle" of Pascal (see Boyer [B3], p. 439). In the arithmetic triangle, 



1 

1 1 

1 2 1 

13 3 1 

1 4 6 4 1 

1 5 10 10 5 1 



each entry (except the first and last in each row) is the sum of the elements to its north and its north-west, 
whereas in the harmonic triangle, 

1/1 

1/2 1/2 

1/3 1/6 1/3 

1/4 1/12 1/12 1/4 

1/5 1/20 1/30 1/20 1/5 

1/6 1/30 1/60 11/60 1/30 1/6 



each entry is the sum of the elements to its south and south-east. The resistance i?„ is the sum of the entries 
in the n-th row: 

n: 1 2 3 4 5 6 7 8 

i?„ : 1 1 5/6 2/3 8/15 13/30 151/420 32/105 

The appearance of the reciprocals of binomial coefficients in (1.1) suggests that we also consider the sum 

E M' (1-2) 

which has the following values: 

n: 1 2 3 4 5 6 7 8 ••• 

Sn-. I 2 5/2 8/3 8/3 13/5 151/60 256/105 83/35 

Of course, these two sequences are linked by the relations 

Rn = -Sn-i and Sn = {n+ l)i?„+i. (1.3) 
n 

In the next section, we shall review some exact results (alternative expressions and generating functions) 
for the numbers i?„ and Sn- In the following section we shall consider asymptotic expansions for these 
numbers. The coefficients in these asymptotic expansions have simple combinatorial interpretations that we 
shall consider in the subsequent section. These combinatorial interpretations will launch us on a tour of 
old and new results in combinatorial enumeration. In the final section, we shall return to the hypercube of 
resistors, and consider the resistance between vertices that are not the endpoints of a long diagonal. 
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2. Alternative Expressions and Generating Functions 



The numbers Rn and Sn have alternative expressions, 



2k 
~k 



(2.1) 



l<fe<n 



and 



Sn — 



n+1 



E 



(2.2) 



0<fc<n 



fc + 1' 



which are equivalent to each other by virtue of (1-3). While having just as many terms as (1.1) and (1-2), 
these sums have simpler summands, and will thus lend themselves more easily to further developments. 



from which (2.2) follows by induction from the base case = 1. The formula (2.1) was given without proof 
by Mullin [M], and was given with an "electrical" proof by Rennie [R], as follows. First, let a current of 

1 ampere flow into a vertex A and out of a vertex B long-diagonally opposite to A. If B is at potential 
0, then A is at potential i?„ volts. Let A' be adjacent to A, and B' long-diagonally opposite to A', and 
therefore adjacent to B. By symmetry, 1/n amperes flows through the 1-ohm resistor from A to A', so A' 
is at potential Rn — 1/n volts. By a similar argument, B' is at potential 1/n volts. Second, reconnect the 
current source so that 1 ampere flows into A' and out of B'. If B' is at potential 0, then A' is at potential 
Rn volts, A is at potential Rn — 1/n volts, and B is at potential 1/n volts. Third, suppose that currents 
of 1 ampere flow in at each of A and A', and out of each of B and B' . By linearity, we may superimpose 
the potentials and subtract 1/n from their sum, putting B and B' at potential 0, and A and A' at potential 
2Rn — 2/n volts. But with this final current distribution, 1 ampere flows from A to B' through the resistors 
of an (n — l)-dimensional hypercube, 1 ampere flows from A' to B through the resistors of another disjoint 
(n — l)-dimensional hypercube, and no current flows through the resistors connecting corresponding vertices 
in the two hypercubes, since by symmetry they at at equal potentials. Thus if B and B' are at potential 0, 
A and A' arc at potential Rn-i volts. We therefore have 2i?„ — 2/n = Rn-i, or i?„ = (l/2)i?„_i -I- 1/n, 
from which (2.1) follows by induction from the base case Rq — 0. Finally, we mention that Sury [S5] proved 
(2.2) by using the integral representation {1 — x)''^~^~^ dx = l/n("^"^) (Euler's beta integral), summing 
the resulting geometric progression inside the integral, and evaluating the resulting integral by a change of 
variable. 

Equations (2.1) and (2.2) allow us to easily derive the generating functions R{z) = X^„>o Rn -2" and 

S{z) = Y,n>o for the sequences Rn and Sn- Indeed, since — log(l— 2) = z+z'^/2+z^/3-\ \-z''/k-\ , 

we see that 2''/k is the coefficient of z'^ in — log(l — 2z). If A{z) = X]n>o -^^ ^^'^ R{^) = X]n>o ^''^ 
the generating functions for the sequences An and respectively, then C{z) = A{Z) B{z) is the generating 
function for the sequence C„ = X]o<fc<n Rn-k, called the "convolution" of the sequences An and As 
a special case, B{z) = 1/(1 — z) is the generating function for the sequence B„ = 1, so that A{z)/{l — z) is the 
generating function for the sequence X^o<fe<rt of partial sums of the sequence An- Thus 2" X^i</j<„ 2'^/A; 
is the coefficient of 2:" in (— log(l — z)) /(I — z/2), so that 



The first proof of (2.2) was given by Staver [S3], who derived the recurrence 



Sn = {{n + l)/2n)Sn-i + l, 



R{z) 




(2.3) 
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Prom (1.3), we see that S{z) = R'{z), so differentiating (2.3) yields 

^"= (l-.)(L./2) + 2(l-W ^°^T^- ^'-'^ 
Generating functions for sums similar to (2.1) and (2.2) have been given by Pla [P]. 



3. Asymptotic Expansions 

The results of the preceding section give exact values of i?„ and S'„ as a rational numbers but they yield 
little insight into the behavior of these sequences for large n. To obtain this insight, we develop asymptotic 
expansions. It will be convenient to use "0-notation" , where 0(/(n)) stands for some function g{n) (possibly 
a different function at each occurrence) such that \g{n)\ < cf{n) for some constant c and all sufficiently 
large n. 

We start with (1.1). Since the binomial coefficients ("]^^) increase as k increases from to [(n — l)/2j, 
then decrease as k increases from \{n — l)/2] to n — 1, the largest terms in (1.1) arc the first and last: 
V("o^) — -'^/("-i) ^ ^- ^^'^ next largest terms are the second and sccond-to-last, which are l/("j^^) = 
I/C-2) = 0{l/n). There are n-4 other terms, and each of these is at most l/i^^^) = I/CZ3) = 0{l/n'^), 
so the sum of all these other terms is also 0(l/n). Thus we have 

Rn = - + 0(^. (3.1) 

This result gives a good estimate of i?„ when n is large. 

We can refine the estimate (3.1) by extracting the second and second-to-last terms, noting that the 
third and third-to-last terms are 0{l/n^), and that each of the remaining n — 6 terms is 0(l/n^), so their 
sum is also 0(l/n^). This yields 

Rn = -(l + +0(^ 

Continuing in this way, we obtain 

^" " n + (^T^ + (n - l)(n - 2) + " ' + (n - l)(n - 2) ■ ■ ■ {n - k) ^ ^ {^)) ^^'^^ 
for any fixed k. 

Equation (3.2) give a sort of asymptotic expansion for i?„, but its content would be clearer if the 
denominator of each term was a power of n, instead of the "falling powers" {n — l)(n — 2) ■ ■ ■ (n — k) that 
appear there. That is, we would like an expansion of the form 

Rn = -(ro + '-^ + ^ + --+'-^+o(^)) (3.3) 

for each A; > 0. It is customary to write 

— ro + — + ^ + -- - + -r + -- - (3.4) 
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as shorthand for the assertion of (3.3) for each A; > 0. The series (3.4) is called an asymptotic expansion; it 
is not convergent for any n, but it allows i?„ to be approximated with an error 0{l/n'^) for any fixed k and 
all sufficiently large n (where the constant hidden in the O- notation depends on k). 

Our task is to determine the coefficients ro,ri,... in (3.4). To do this we shall expand each term 
k\/{n— l)(n — 2) • • • (n — fc) in (3.2) into a series of negative powers of n, 

k\ yr^kltk 



(n-l)(n-2)-..(„-t) 

then sum the contributions to r; for each k < I. First we need to find the numbers tk,i in the expansion (3.5). 
These arc what have come to be called the "Stirling numbers of the second kind" , for which we shall use the 
notation suggested by Knuth [K, p. 65]: tk,i ~ {i}- These numbers were introduced by James Stirling in 
the Introduction to his Methodus Differentialis [S4] in 1730. He defined them as the numbers that expand 
a power 0' of 2; as a linear combination of the polynomials z, z{z — 1), ... , z{z — 1) ■ ■ ■ {z — I + 1): 



z' 



PU(2-i)...(2-fc+i), 



Q<k<l 



and he gave a table for 1 < fc < ? < 9. The number | ^ | has a simple combinatorial interpretation: it is the 
number of ways to partition the I elements of the set L = {1, . . . , / } into k blocks (non-empty subsets of L 
that are pairwise disjoint and whose imion is L). For I = 3, for example, we have one partition {{1,2,3}} 
into one block, three partitions {{1}, {2, 3}}, {{1, 2}, {3}} and {{1, 3}, {2}} into two blocks and one partition 
{{1},{2}, {3}} into three blocks; thus {^}=1, {^}=3 and {3} = 1. At the end of the introduction, 
Stirling gives the expansion 



(2 + l)(2 + 2)...(2 + fc) ' \k\ z^' 



which, upon substitution of — n for z, gives (3.5) in the form 



kx ^^,rn 1 



{n-\)(n-2)---(n-k) ^ ' \k] n^' 

These expansions are in fact convergent for fixed A; > 1 and for |2;| < 1/fc or n > fc (though Stirling did not 
distinguish convergent expansions, such as these, and asymptotic expansions, such as (3.4)). Applying this 
result to each term in (3.2) gives the desired asymptotic expansion: 

1>Q \0</s<( 



Thus the coefficients ri we sought are given by 



0<fe<i 



(It will doubtless have occurred to the reader that if there are "Stirling numbers of the second kind" , 
there should also be "Stirling numbers of the first kind" . Indeed there are, and they were also introduced 



5 



by Stirling [S4]. He defined them as the numbers that expand z{z + 1) • • • + Z — 1) as a Unear combination 

of the polynomials z, z'^, . . . ,zK Nowadays it is more common to define them as the absolute values of the 
numbers that expand z(z — 1) ■ ■ ■ [z — I + 1) as a linear combination of the polynomials z, z"^, . . . , z'; in the 
notation of Knuth [K, p. 65]: 



0<l<l 

Stirling again gave a table for 1 < < Z < 9, and the expansion 

"n 1 

k -^-^\■■ 

l>k 



z\ 



z{z+l)---{z + k-l) 



These numbers too have a simple combinatorial interpretation: ^ is the number of permutations of I 
elements that have k cycles. For I = 3, for example, we have two permutations (123) and (132) with one 
cycle, three permutations (1)(23), (12)(3) and (13)(2) with two cycles and one permutation (1)(2)(3) with 
three cycles; thus [\] = 2, [\] = 3 and [3] = 1.) 

We can find a similar asymptotic expansion for Sn- Again noting that the largest terms in the sum (1.2) 
are the first and the last, we obtain 

5„ = 2 + O 



Generalizing this as before yields 



n n{n — 1) 



+ 



k\ 



n{n — 1) • • • (n — fc + 1) 



O 



,fe+i 



Applying (3.5) to each term and summing the contributions for each negative power of n, we obtain 

Sn Si S2 , Sfe 

— ~ So + — + — + ••• + — + ■ 

2^ vt ■n-^ •n"' 



n 



(3.7) 



where 



si 



E 

0<k<l 



(A: + l)!. 



The coefficients ri and si have simple combinatorial interpretations that we shall study in the following 
section. 



4. Preferential Arrangements 

In this section we shall study the numbers r/ and s;. Our model for this study will be a collection of 
results concerning the "exponential numbers" di, given by 



*= E III- (") 



0<fe<( 



These numbers have a simple combinatorial interpretation: di is the number of ways to partition the set 
{!,...,/} into any number of blocks. For I = 3, for example, we have seen that there is one partition into 
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one block, three partitions into two blocks and one partition into three blocks; thus ^3 = 1 + 3 + 1 = 5. We 
have the table 

Z:01234 5 6 7 8 
d; : 1 1 2 5 15 52 203 877 4140 • • • . 

(The sequence di is AOOOllO in Sloan's On-Line Handbook of integer Sequences [S2].) 

There are three aspects of the exponential numbers that are of particular interest to us: a recurrence, 
a generating function and an expression as the sum of an infinite series. The recurrence is 

di=Si+ Cl^)'i>'^ (4.2) 



0<k<l-l 



where Si is 1 for Z = and for all other values of I. This recurrence allows di to be computed from the 
previous values do,di, . . . ,di-i. 

The "exponential generating function" d{z) = J2i>o -^V^' given by 

d{z) = e«''-\ (4.3) 

where the term "exponential" refers to the factor in the defining sum (in contrast to the "ordinary" 
generating functions that we used in Section 2). 

The expression as an infinite sum is 

dt=-j:^, (4.4) 

e n! 

n>0 

where e = 2.7182 ... is the base of natural logarithms. Since di is expressed as a finite sum in (4.1), it may 
not be clear what advantage there is to (4.4). But the finite sum involves the Stirling numbers of the second 
kind, whereas the infinite sum involves only powers and factorials. Furthermore, the sum (4.4) has a simple 
probabilistic interpretation: di is the l-th moment Ex[A'''] of a Poisson-distributed random variable N with 
mean A = 1 (since Pr[A'' = n] = e~^ A"/n! = 1/en! for such a random variable). 

The exponential numbers were mentioned in 1934 by Boll [Bl, B2], and arc on that account sometimes 
called the "Bell numbers". But (4.2), (4.3) and (4.4) were all given earlier: (4.2) in 1933 by Touchard [T] 
(in an equivalent "umbral" form), (4.3) in 1886 by Whitworth [W, Proposition XXIV, p. 95] and (4.4) in 
1877 by Dobinski [D] (who actually only gave the cases 1 < / < 8; but it is clear from his derivations that di 
satisfies the recurrence (4.2)). 

In his marvelous book Asymptotic Methods in Analysis, N. G. de Bruijn [B5, Section 3.3] derives the 
asymptotic expansion 

1 , ■ ^ do di dk 

— > fc!~l + — + — + ••• + -^ + ---, 

a<k<n 

and then says that it is "only for the sake of curiosity" that he mentions that the coefficients d;, given by 
(4.1), have a combinatorial interpretation. One of our goals in this paper is to pursue this curiosity; our 
motto is: whenever the coefficients in an expansion are integers, look for a combinatorial interpretation! 

The coefficients r;, given by (3.6), also have a simple combinatorial interpretation: they are the number 
of ways of ranking / candidates, with ties allowed; that is, the I candidates are first to be partitioned 
into equivalence classes, then the equivalence classes are to be linearly ordered. This interpretation follows 
from those of |^| and k\, where k is the number of equivalence classes in the partition. Because of this 
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interpretation, rj is called the number of preferential arrangements of / elements. For I = 3, for example, 

the one partition into one block can have its block ordered in one way, each of the three partitions into two 
blocks can have its blocks ordered in two ways and the partition into three blocks can have its blocks ordered 
in six ways; thus r3 = l- l + 3- 2 + l- 6 = 13. Wc have the table 

;:01234 5 6 7 8 

n : 1 1 3 13 75 541 4683 47293 545835 • • • . 

(The sequence r; is A000670 in Sloan [S2].) We shall derive the recurrence 

ri=di+ Yl Uj^fc' (4-5) 

0<k<l-l ^ ^ 

the exponential generating function (defined by r{z) = X^oqH z'//!) 

K^) = ^ (4.6) 

and the summation expression 

= 2 ^ ^- (4-7) 

n>0 

We begin by deriving the recurrence (4.5). For Z > 1, we can construct a preferential arrangement on 

I candidates by first choosing the number k of candidates tied in the top equivalence class (with k in the 
range 1 < k < I), then choosing (in one of (^) ways) the candidates in this class, and finally choosing (in 
one of Ti-k ways) a preferential arrangement of the remaining I — k candidates. This gives the recurrence 

l<k<l ^ ' 



E 

0<j<(-l 



E 

0<fc<i-l 

E 



l-k 



n-j-i 



rk, 



0<k<l-l 



rk, 



where we first made the substitution k = j + 1, then the substitution j = I — 1 — k, and finally used the 
identity (jij.) = {j.)- This equation holds for I > 1; since ro = 1, we obtain (4.5) for I > 0. 

Next we shall derive the exponential generating function (4.6). Adding r; to both sides of (4.5) yields 

2ro = 6i+ ([) ^J- (4-8) 

o<j<i ^-'^ 

Multiplying both sides of this equation by z''/l\ and summing over Z > 0, we obtain 



M^) = i+E^ E G 



^>0 0<j<l 



^ j! fc! 

j>o fe>0 

= l + e^r(^), (4.9) 
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where we have made the substitution k = l—j and used the identity = X^fe>o z'^/kl. Solving this equation 
for r{z) yields (4.6). 

Finally we shall derive the summation expression (4.7). To do this, we rewrite the exponential generating 
function r{z) from (4.6): 

^riz^ _1 1 




Since the coefficient of z^/l\ must be the same on both sides of this equation, we obtain (4.7). 

The name "preferential arrangement" was introduced by Gross [G], as was the summation expression 
(4.7). The numbers r/ (with a different combinatorial interpretation involving trees), the recurrence (4.5) 
and the generating function (4.6) were given by Cayley [CI] in 1859; the combinatorial interpretation we 
have used is implicit in 1866 by Whitworth [W, Proposition XXII, p. 93] (Whitworth shows that the terms 
I ^ I fc! for fixed k have the exponential generating function (e^ — 1)'^; summation over A; > then yields 
1/(1 -(e^-l))= 1/(2 -e^).) 

We turn now to the numbers s;, which also have a simple combinatorial interpretation: s; is the number 
of ways of ranking I candidates, with ties allowed, and with a "bar" that may be placed above all the 
candidates, between two equivalence classes of tied candidates, or below all the candidates. Thus we may 
call s; the number of barred preferential arrangements of I elements. If there are k equivalence classes of tied 
candidates, there arc fc + 1 positions for the bar. For I = 3, for example, the one preferential arrangement 
with one block has two positions for the bar, each of the six preferential arrangements with two blocks has 
three positions for the bar, and each of the six preferential arrangements with three blocks has four positions 
for the bar; thus S3 = 1 • 2 + 6 • 3 + 6 ■ 4 = 44. We have the table 

^:01234 5 6 7 8 

s; : 1 2 8 44 308 2612 25988 296564 3816548 ■ • • . 

(The sequence s; is A005649 in Sloan [S2].) 

In lieu of a recurrence for the numbers s„, we shall derive a formula expressing them in terms of the 
numbers rf. 

0<k<l ^'^^ 

We shall also derive the exponential generating function (defined by s{z) = J2i>o -^'/^O 

= (4-11) 
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and the summation expression 



n>0 



For Z > 0, we can construct a barred preferential arrangement on / candidates by first choosing the 
number k of candidates above the bar (with k in the range < fc < then choosing (in one of (^) ways) the 
candidates above the bar, then choosing (in one of ways) a preferential arrangement of these candidates, 
and finally choosing (in one of ri-k ways) a preferential arrangement of the remaining I — k candidates. This 
gives the formula (4.10). Next, multiplying both sides of (4.10) by z^/l\ and summing over I > yields 



;'>n n<-ii.<-; ^ 



l>0 0<k<l 



I k \ " ' j 

k\ ^ j\ 

k>0 j>0 ■' 

= r{zf, (4.13) 

where we have made the substitution k = I — j. Substituting (4.6) in this equation yields (4.11). Finally, 
reasoning similar to that used to derive (4.7) leads to (4.12). 

Before concluding this section, let us derive two more identities relating r; and sf. 



0<k<l 



and 

ri+ri+i=2si. (4.15) 

These can be given direct combinatorial proofs (and the reader may enjoy finding these), but we shall use two 
different methods that are often useful when dealing with sequences that have explicit exponential generating 

functions. 

To prove (4.14), we use the notion of "binomial convolution". Suppose that a(z) — X/;>o'^'^V^' s-nd 
b{z) = X]i>o /l^- are the exponential generating functions for the sequences ai and hi, respectively. Then 

«WK^) = E^E^ 

fe>o ■ j>0 •'' 



l>0 0<k<l 

where we have made the substitution j = l — k. Thus c{z) = a{z) b{z) is the exponential generating function 
for the sequence c; = X]o<j<( (') '^j ^i-j^ which is called the binomial convolution of the sequences a; and 
hi, and denoted (a * b)i. (We have already encountered binomial convolutions twice in this section: once to 
derive (4.9) from (4.8), where the convolution can be expressed as s = S + r * v (the sequence vi = 1 for 
all / > 0, and has exponential generating function v{z) = e^), and again to derive (4.13) from (4.10).) We 
shall also need the fact that a'{z) (where the prime indicates differentiation) is the exponential generating 
function for the sequence a;+i, which we shall denote a^. 
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To derive (4.14), we may now observe that r'{z) = e^/(2 — e^)^ = s{z). Thus n+i = r'f = {s * v)i, 
which yields (4.14). 

To derive (4.15), we note that r'{z) = e^/(2 — e^)^ impUes that r{z) satisfies the differential equation 

r'{z) + r{z) = 2r{zf. (4.16) 

(This differential equation, together with the initial condition r(0) = 1, uniquely determines r(z). In fact, it 
is an example of a "Riccati equation", which can be solved to analytically for r{z).) Substituting (4.13) in 
(4.16), wc obtain (4.15). (We note that (4.15) can also be obtained from (4.7) and (4.12).) 

That the numbers si, defined by (3.7), have the exponential generating function given in (4.11) was given 
as an exercise (without proof or reference) by Comtet [C2, p. 294, Ex. 15]. Our combinatorial interpretation 
of these numbers in terms of barred preferential arrangements seems to be new. 



5. More Resistances 

We mentioned in the introduction that Singmaster [SI] posed in 1978 the problem of determining i?„. 
What we did not mention then is that he asked not only for i?„, but for Rn,k- the resistance between two 
vertices at distance k (for 1 < fc < n) in an n-dimcnsional hypcrcube of 1-ohm resistors. Ronnie's solution 
[R] to Singmaster's problem covered (by various arguments) the cases k = 1,2 and k = n,n — l,n — 2. In 
this section we shall give (by a single argument) a complete solution to Singmaster's problem: for < A: < n, 

2 11^ fn\ , , 

^n.. = - E T^T^F E (J- (5-1) 

0<j<fe-l \ 3 I j+l<i<n ^ ^ 

Our solution, like that of Rennie, is based on the principles of symmetry and superposition. 

Consider the situation in which a current of 1 ampere flows out of a vertex A, while currents of 1/(2" — 1) 

amperes flow into each of the 2" — 1 other vertices. Symmetry ensures that all (") vertices at distance j from 
A are at the same potential. Call this potential Uj volts, where C/q = 0. There are l/n("T^) 1-ohm resistors 
connecting vertices at potential Uj to vertices at potential ?7j+i, and a total current of ^^.j.^^^^^^ (") /(2" — 1) 
amperes flows through them. By Ohm's law, 

'n 



n{^-^\ 2"-l ^ \i 

V 7 / 7 + l<i<n ^ 



3 I j+l<i<n 

and thus 



1 1 I _ 



n ^ ("7^) 2" - 1 ^ V ^ 

0<j<fe-l \ ] I j+l<i<n 



Now let B be a vertex at distance k from A, and consider the situation in which a current of 1 ampere flows 
into B and currents of 1/(2" — 1) amperes flow out of each of the 2" — 1 other vertices. In this situation there 
is again a potential difference of IJk volts between A and B. By superposition, if a current of 1 + 1/(2" — 1) 
amperes flows into B and out of A, there will be a potential difference of 2(7^ between these vertices. Again 
using Ohm's law, we have 



r. 2 1 ^ _l \_ ^ (n 

l + l/(2"-l)n ^ f"-i) 2"-l ^ 

0<j<k-l V J / j+l<i<n 
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which yields (5.1). 

Let us first check that (5.1) agrees with (1.1) for k = n. We have 



0<j<n-l V J / j+l<j<rt ^ 
1 1 1 ^ /n 



n ^ ("7^) 2" ^ v« 

0<j<n-l \ j J j+l<i<n 

n 



1 ^ ^ 

' 0<j<n-l y j ' j + l<i<n ^ 

1 >r- 1_ sr- 



0<i<n-l V J / 0<i<i 



= - V V 

n ^ 2" ^ 

= - V ^ 

n ^ ("7^) 

0<j<n-l V J / 



where we have used the identities ("^ ) = and (") = („" j). 

Next, let us show that Rn,k is an increasing function of k (that is, that \7kRn,k — Rn,k — Rn,k-i > 
for 1 < A; < n). From (5.1) we have 

Vfe^in,. = t;!^ E (•)' (5-2) 

and the expression on the right-hand side is obviously positive. 

Finally, let shall show that Rn,k is a concave function of k (that is, that V^Rn,k = Rn,k — '^Rn,k-i + 
Rn,k-2 < for 2 < A: < n). From (5.2) we have 



\ \k-\) k<i<n ^ ' \k-2) k-l<i<n 



The expression on the right-hand side is obviously negative for k < n/2 + 1, since in this case we have 
(^2) < (r!) and Ek<i<n G) < T.k-l<^<n (") ' ^"^^ ^ase k>n/2 + 1, we factor (^Z^) out of the 
expression in parentheses, then move all but the first term from the second sum into the first sum: 





( k-l 


„2n-l (-1) 1 


\n-k+l 




1 2k-n-2 


^2-1 (-1) 1 


1 n-k+1 



k<i<n k—l<i<n 



k<i<n 



n\ In 
i)~\k-l 
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Thus what remains to be proved is that the expression in parenthesis is negative; that is, that 

2fc — n — 2 \ ^ fri\ f n 



- y 
1 ^ 



n — k + l , \i J \k — l 



< 



for k > n/2 + 1, or equivalently, by the substitution k = n — j and the identity (") = („" J, 

n — 2 j — 2 ■r-^ f n\ f n 



7 + 1 ^ \i J V J + 1 

for < j < n/2 - 1. 

To prove (5.3), we observe that < i < n/2 — 1 impHes that 

n\ ^ z + 1 / n 



n — i\i + l 

Since (i + l)/(n — i) is an increasing function of i, we have 

\i)-\n-j) + 
for i < j < n/2 — 1. Thus we may bound the sum in (5.3) by the sum of an infinite geometric series, 

0<i<j ^ ^ 0<i<j 

J + 1 



■ I -I \ m / 

J + l\ ( n 



2 - 2j - 1 \j + 1, 

This inequahty proves (5.3), and thus completes the proof that Rn,k is concave. 



6. References 

[Bl] E. T. Bell, "Exponential Polynomials" , Annals of Mathematics, 35:2 (1934) 258-277. 

[B2] E. T. Bell, 'Exponential Numbers", American Mathematical Monthly, 41:7 (1934) 411-419. 

[B3] C. B. Boyer, A History of Mathematics, John Wiley & Sons, 1968. 

[B4] E. E. Brooks and A. W. Poyser, Electricity and Magnetism: A Manual for Advanced Classes, Longmans, 
Green, 1914. 

[B5] N. G. de Bruijn, Asymptotic Methods in Analysis, North-Holland, 1958 (reprinted by Dover, 1981). 
[CI] A. Cayley, "On the Analytical Forms Called Trees— Part 11" , Philosophical Magazine, 18 (1859) 374-378. 
[C2] L. Comtet, Advanced Combinatorics, D. Reidel Publishing, 1974. 

[D] G. Dobinsky, "Summierung der Reihe ^n^/n! fiir m = 1,2,3,4,5,...", Archiv fiir Mathematik und 
Physik (= GruneH's Archiv), 61 (1877) 333-336. 

[G] O. A. Gross, "Preferential Arrangements", American Mathematical Monthly, 69:1 (1962) 2-4. 

13 



[I] E. L. Incc, Ordinary Differential Equations^ Dover, 1956. 

[J] A. A. Jagers, "Solution to Elementary Problem E 2620" (also solved by T. Morley and D. Zave), 

American Mathematical Monthly, 85:2 (1978) 117-118. 
[K] D. E. Knuth, The Art of Computer Programming, v. 1: Fundamental Algorithms, Addison- Wesley, 1968. 
[M] A. MuUin and D. Zave, "Elementary Problem E 2620", American Mathematical Monthly, 83 (1976) 740. 

[P] J. Pla, "The Sum of the Inverses of the Binomial Coefhcients Revisited", Fibonacci Quarterly, 35: 4 
(1997) 342-345. 

[R] B. C. Rennie, "Solution to Problem 79-16", SIAM Review, 22 (1980) 504-508. 
[SI] D. Singmaster, "Problem 79-16", SIAM Review, 21:4 (1978) 559. 

[S2] N. J. A. Sloane (Ed.), The On-Line Encyclopedia of Integer Sequences, published electronically at 
http : //www. research. att . com/ ~nj as/sequences , 

[S3] T. B. Staver, "Om summasjon av potenser av binomiaalkoeffisientene" , Norsk Matematisk Tidsskrift, 
29 (1947) 97-103. 

[S4] J. Stirling, Methodus Differentialis, 1730 (Translated and annotated by I. Tweddle, Springer, 2003). 

[S5] B. Sury, "Sum of the Reciprocals of the Binomial Coefficients", European J. Combinatorics, 14:4 (1993) 
351-353. 

[T] J. Touchard, "Proprietes arithmetiques de certains nombres recurrents", Ann. Soc. Sci. Bruxelles, 53 
(1933) 21-31. 

[W] W. A. Whitworth, Choice and Chance, Deighton, Bell and Co., 1886. 



14 



